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Abstract 

Analyticity of gluon and Faddeev-Popov ghost propagators and their form 
factors on the complex momentum-squared plane is exploited to continue ana- 
lytically the ultraviolet asymptotic form calculable by perturbation theory into 
the infrared non-perturbative solution. We require the non-perturbative multi- 
plicative renormalizability to write down the renormalization group equation. 
These requirements enable one to settle the value of the exponent characteriz- 
ing the infrared asymptotic solution with power behavior which was originally 
predicted by Gribov and has recently been found as approximate solutions of 
the coupled truncated Schwinger-Dyson equations. For this purpose, we have 
obtained all the possible superconvergence relations for the propagators and 
form factors in both the generalized Lorentz gauge and the modified Maximal 
Abelian gauge. We show that the transverse gluon propagators are suppressed 
in the infrared region to be of the massive type irrespective of the gauge pa- 
rameter, in agreement with the recent result of numerical simulations on a 
lattice. However, this method alone is not sufficient to specify some of the 
ghost propagators which play the crucial role in color confinement. Combin- 
ing the above result with the renormalization group equation again, we find 
an infrared enhanced asymptotic solution for the ghost propagator. The cou- 
pled solutions fulfill the color confinement criterion due to Kugo and Ojima 
and also Nishijima, at least, in the Lorentz-Landau gauge. We also point out 
that the solution in compatible with color confinement leads to the existence of 
the infrared fixed point in pure Yang-Mills theory without dynamical quarks. 
Finally, the Maximal Abelian gauge is also examined in connection with quark 
confinement. 
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1 Introduction 



Quantum chromodynamics (QCD), a non-Abelian gauge theory P with a gauge group 
SU{3), is now regarded as the fundamental theory for describing strong interactions 
among quarks, gluons and hadrons as their bound states. Since asymptotic freedom [2j 
exists in QCD with Nf quark flavors less than 33/2, one can predict the ultraviolet 
(UV) behavior of QCD by perturbation theory in a reliable way. In order to discuss the 
infrared (IR) behavior, on the other hand, we must solve the strong coupling problem 
beyond the perturbation theory. In the IR region, indeed, the perturbation theory 
itself cannot be a reliable method of calculations. Hence, we need inevitably non- 
perturbative methods. However, this does not necessarily mean that the perturbation 
theory cannot provide any information on the IR behavior of QCD, as will be shown 
in this article. 

One of the most direct methods to study the non-perturbative features of quantum 
field theory would be to write down and solve the Schwinger-Dyson (SD) equation. 
The SD equation is a set of self-consistent equations among Green functions and 
constitutes a series of exact equations. However, the SD equations consist of the 
infinite number of simultaneous non-linear integral equations. In actually solving the 
SD equation, therefore, one needs to perform a certain truncation which can not be 
justified within the framework. Nevertheless, the approximate solutions have become 
an efficient tool for studying non-perturbative problems which cannot be reached by 
perturbative method, e.g. chiral symmetry breaking. See e.g., reviews [SlIlllS]- 

So far, the extensive studies of the SD equations in QCD have almost been re- 
stricted to the SD equation for quark propagator (referred to as quark SD equation 
hereafter), see Figure ^a). In order to study gluodynamics,^ we must solve the cou- 
pled SD equations for the gluon and Faddeev-Popov (FP) [7] ghost propagators, see 
Figure ^b),(c). Of course, solving them is more difficult than solving the quark equa- 
tion alone. This is mainly because the SD equations for gluon and ghost are truly 
non-linear coupled equations. This is in sharp contrast with the quark SD equation 
studied so far. 

An origin of this non-linearity comes from the treatment in which both propagators 
are incorporated on equal footing, see Figure ^b),(c). In the conventional treatment 
of the quark SD equation, the gluon propagator was replaced by hand with the free 
propagator (ladder approximation) or dressed propagator in the one-loop resumed 
level (improved ladder approximation) and consequently the problem was reduced to 
solving the quark propagator alone. ^ Another origin of non-linearity comes from the 
fact that gluon and ghost are bosons. The quark SD equation can be decomposed into 
two equations for the mass function and the wave-function renormalization function. 
The equation for the quark mass function can be linearized as a good approximation 
or in the rigorous sense (bifurcation technique) due to a fact that quark is a fermion. 
The general solution for the mass function can be obtained by adopting either method 
of linearizations. Once the general solution is obtained, it can be used to study the 

^In this article, we deal only with the manifestly Lorcntz covariant formulation of Yang-Mills 
theory. Therefore, we do not discuss the non-covariant gauge, such as the Coulomb gauge and axial 
gauge ini- We do not discuss the Fock-Schwinger gauge in this article either. 

^For the four-dimensional QED, the coupled SD equations for electron propagator and photon 
propagator in the Landau gauge were solved a decade ago under the bare vertex approximation in 
an article |H] and subsequently by Bloch and Pennington PI in a numerical way. 
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whole range of momenta. This fact facihtates considerably the analysis of the quark 
mass function. 




Figure 1: The coupled Schwinger-Dyson equations for three propagators in QCD with 
the conventional Lorentz gauge fixing, (a) quark equation, (b) ghost equation, (c) 
gluon equation. Here S denotes the full quark propagator, A the full ghost propagator 
and D the full gluon propagator, while F denote the four types of vertices. The pure 
Yang-Mills part is enclosed by the broken line. 

In recent several years, extensive analytical studies [TUllTTllT^IT^inilT^IT^ITTl 
UHl 11111201 1211 1221 123 of the coupled SD equations for gluon and ghost propagators 
were performed to obtain the infrared (IR) and ultraviolet (UV) asymptotic solu- 
tions in the SU{Nc) Yang-Mills theory restricted to the Landau gauge. The lead- 
ing term of the UV asymptotic solution recovers the one-loop resumed perturba- 
tion theory due to asymptotic freedom. In fact, the solutions for the Euclidean 
gluon propagator D^^{Q'^) = S^^{6^^ - Q^,Qu/Q'^)Dt{Q'^) and ghost propagator 
A^f (g2) = 5"^^ /\fp{Q'^) (A, 5 = 1, ■ ■ ■ , A/-^ - 1 for SU{N^)) have the expected UV 
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forms 

^T(g') = ^ (log QT, A^p(Q2) = ^(logQY, (1.1) 

where C and D are Q^-independent constants and the exponents 7 and 5 are inde- 
pendent of the number of colors N^. in SU{Nc) Yang-Mills theory at one-loop level, 
i.e., 7 = -1| and (5 = 

A remarkable characteristic of gluodynamics claimed by von Smekal, Hauck and 
Alkofer ^0] and subsequently confirmed by Atkinson and Bloch in a simpler 
framework is that the Euclidean gluon propagator behaves like a power of the mo- 
menta in the IR limit of vanishing momenta — ^ 0: 

Dt{Q') = A(Q2)2-i^ (1.2) 

and the gluon propagator is suppressed in the IR regime for k > 0, while the Euclidean 
ghost propagator is enhanced to be more singular than the free one in the IR limit 

AMQ') = B{Q^r^-\ (1.3) 

where A and B are Q^-independent constants. This result is surprisingly due to 
dominant contribution from the ghost loop in IR region. However, the value of the 
exponent k varies in the range < k < 1 depending on the way of truncation 
[niliailHlliaiiaiiainillHlliniEDliniEai^ At present, the precise value of k is not 
known. Therefore, it is still under the debate whether the gluon propagator converges 
to a finite non-zero value (k = 1/2) or vanishes (k > 1/2) in the IR limit — 0.^ 
It should be remarked that such behaviors of gluon and ghost propagators were first 
predicted by Gribov |23j. His result corresponds to k = 1: 



DTm-T^r^^ ^FP{Q'') (1.4) 



where the mass scale Mq is called the Gribov scale. 

A serious problem is that at present one does not have any analytical method for 
connecting in the intermediate region two asymptotic solutions to obtain the complete 
solution over the whole range of momenta. Only the numerical studies have succeeded 
in obtaining the complete solution in this sense. Therefore, there is no guarantee that 
a specific IR asymptotic solution obtained in the analytical study can be a piece of 
the complete solution which must reduce to the UV asymptotic solution in the large 
momenta. There is also a problem of uniqueness of the solution jlOj, since the IR 
asymptotic solution is obtained based on a specific Ansatz for the solution. Therefore, 
one cannot judge if a set of IR power solutions p.2|) and p.3|l is unique for the coupled 
SD equations. The other problem is that all the SD equations aforementioned have 
been solved in Euclidean space. Indeed, the obtained Euclidean solutions are very 
helpful for comparison with the results of numerical simulations on a lattice. However, 
one needs to perform analytic continuation in order to go back to the Minkowski 



■^On the lattice, there is a general argument [21] which supports the vanishing of the gluon 
propagator at Q = 0. However, this is not necessarily the case in the continuum limit due to the 



renormalization effect, see pp. 146-147 of )25| . 
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spacetime. The analytic continuation to the Minkowski region is not necessarily 
ensured for the Euclidean solutions with the asymptotic behaviors ()1.2|) and ()1.3|) . 

The main purpose of this article is to study the relationship between quark, gluon 
and ghost propagators in Euclidean region and the spectrum of QCD in the Minkowski 
region in connection with color confinement. In this article, we show in Yang-Mills 
theory with the gauge group SU (Nc) that the analyticity of the full propagator on the 
complex momentum plane plays the role of connecting the perturbative UV asymp- 
totic solution into the non-perturbative IR asymptotic solution and vice versa. This 
study resolves the question as to what the IR asymptotic solution with power behavior 
(characterized by an exponent k) means for the spectrum with respect to the funda- 
mental fields. We suggests which type of IR asymptotic solutions should be searched 
for in the SD equation, while the UV asymptotic solution can be reliably predicted 
by perturbation theory due to asymptotic freedom. In other words, the precise value 
of hi will be predicted, if the solution behaves like a power of the momenta. 

In this article, we examine two types of manifestly Lorentz covariant gauges: first 
the Lorentz gauge in the most general form pT) | l27 | l28 ^ l29j , so called the Curci-Ferrari 
gauge, and then the modified version j22ll22lElEiESlEZll2niEillillIlliail3lll 
of the Maximal Abelian (MA) gauge [2011^ • In the conventional Lorentz gauge, 
all the approximate solutions of the coupled SD equations for the gluon and ghost 
propagators have been obtained only in the Landau gauge. The SD equation approach 
in the MA gauge in the same level as the Lorentz gauge is still missing in the MA 
gauge, but see jUj for a preliminary result. 

This article is organized as follows. In section 2, we introduce the generalized 
Lorentz gauge and the modified MA gauge for later convenience. For the above 
purpose, we analyze the propagators in the complex k"^ (momentum squared) plane 
in section 3. We require non-perturbative multiplicative renormalizability for writing 
down the renormalization group (RG) equation for the propagators in the complex 
k"^ plane. The analyticity of the propagators except for the positive real axis and 
the asymptotic behavior for large enables us to specify the IR asymptotic form 
(|P| 0) such that it can be analytically continued to the UV asymptotic form 
{\k'^\ — > oo) calculable by perturbation theory. 

In this article, we investigate which type of IR asymptotic solution of the SD 
equations is consistent with the general principles of quantized gauge field theories: 

1. Non-perturbative multiplicative renormalizability 

2. Analyticity 

3. Spectral condition (and completeness) 

4. Poincare group structure 

without relying on the approximate solutions of the SD equation. 

The former half of this article is deeply indebted to an idea of old articles by 
Oehme and Zimmerman [^ITF] and subsequent articles by Oehme and his collabora- 
tor |lHlliniEniEIlE2llMlEllEniEniEZl • A concept of the superconvergence is crucially 
important also in this article. The superconvergence relation is nothing but a sum 
rule for the spectral function of the fundamental field propagator which is defined 
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in the Minkowski region.^ We reconsider anew the relationship between the super- 
convergence relation and the Euclidean propagator, which is to be compared with 
the solution of the coupled SD equations, paying special attention to the IR asymp- 
totic solution. We explicitly obtain all the possible superconvergence relations in the 
generalized Lorentz gauge in section 4 and the modified MA gauge in section 7. In 
the generalized Lorentz gauge, we show that the superconvergence relation holds for 
gluon. Consequently, the Euclidean gluon propagator has the power-series expansion 
in and goes to a Q"^ -independent constant, i.e., Dt{Q^) ~ const, in the IR limit 
— * 0. Therefore, this result suggests that the IR exponent k, characterizing the 
Euclidean solution must take the value 



irrespective of the value of gauge fixing parameter. 

A question is whether or not the IR asymptotic solution for the ghost propagator 
exhibits power behavior, just as the recent approximate solution of the truncated SD 
equation suggests. So far, one has payed little attention to the superconvergence for 
ghost. We claim that the superconvergence does not hold for ghost in the pure Yang- 
Mills theory (without dynamical fermions) against the previous result jnEj- In fact, 
we elucidate that non-existence of the superconvergence for ghost is consistent with 
realizing color confinement in QCD. At the same time, we point out the subtleties of 
using the Landau gauge in the Lorentz type gauge fixing, although all the approximate 
solutions of the coupled SD equations for the gluon and ghost propagators have been 
obtained only in the Landau gauge. In this article we shall provide a unified view 
among Euclidean solutions of the SD equations, color confinement and the spectrum, 
especially the mass gap of Yang-Mills theory and QCD. 

In section 5, the IR asymptotic form in Euclidean region {k"^ < 0: negative real 
axis) is related to the spectral property in the Minkowski region {k"^ > 0: positive 
real axis) by way of the spectral representation which follows from the spectral condi- 
tion, Poincare covariance and completeness of the intermediate states. We show that 
Dj-{0) < oo is consistent with the existence of a massive pole in the gluon propagator 
in the Minkowski region.^ Moreover, we show that, for gluon with massive spectrum, 
the gluon propagator and the gluon form factor have the well-defined power-series ex- 
pansions in for small Euclidean momenta Q^, suggesting the uniqueness of the IR 
asymptotic solution with power behavior. In the Landau gauge, this result does not 
contradicts with the existence of the gluon condensate (Ajj) of mass dimension two 

which has been extensively studied by several groups inainiiinsiiM!, see ^ for more 

details on the physical meaning of dimension two condensate. The gluon condensate 
(A^) is a special limit of the ghost-gluon condensate of mass dimension 2 which is 
on-shell BRST invariant as proposed in [B31inn] and calculated recently in [SSj . 

In the latter half of this article, we shall focus on the following issues. For the FP 
ghost propagator, however, the analyticity approach alone is not sufficient to provide 

^The various sum rules were devised in QCD to investigate the hadronic properties, see [58I59IS0] 
for reviews. The sum rule in this article is restricted to that for the fundamental field, rather than 
the composite hadronic field operators. 

^It may happen that this constant DxiO) vanishes. This case implies k = 1, the Gribov limit. 
This possibility can not be excluded from the general considerations only. However, such a case is 
realized only when the subtle cancellation occurs between two terms. The argument in more detail 
will be given later. 
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precise form of the IR asymptotic solution. This is reasonable for our result to be 
consistent with the result of SD equation, since the superconvergence for ghost in the 
Landau gauge contradicts with the solution of the SD equations mentioned above, 
i.e., enhancement in the IR region. In section 6, therefore, we make use of the RG 
equations for the gluon and ghost propagators in the IR region. First, we point 
out that the RG equation for the ghost propagator has the IR asymptotic solution 
with power behavior, if the RG functions in the RG equations have the power series 
expansion around the IR limit of the coupling constant g'^{fi = 0), see Figure Here 
we use the fact that the gluon propagator has the IR asymptotic solution with power 
behavior, which follows from the superconvergence argument above, and that the 
P function is common to both RG equations. Thus the power solution Dt{Q^) = 
A{Q'^)°'~^ for the gluon propagator yields the power solution App^Q^) = B{Q'^Y~^ 
for the ghost propagator whose IR critical exponent is characterized by the anomalous 
dimension of the ghost. 

Next, we must answer how to obtain the exponent (3 in the IR asymptotic solutions 
for the ghost propagator. The non-renormalization of the renormalization factor Zi 
for the gluon-ghost-antighost vertex in the Landau gauge yields the relation (3 = 
—a/2 = —K, by supposing the existence of an IR fixed point, i.e., a finite value 
of as{0) := suggested from the solution of the SD equation. Thus the ghost 

propagator in the Landau gauge must behave like 

Appm = m-'^' (1.6) 

in the IR region, provided that k, = 1/2. This result shows that the color confinement 
criterion due to Kugo and Ojima 

JimJQ2A^p(g2)]-i ^ I + ^^0) = (1.7) 

is fulfilled in the Landau gauge Yang-Mills theory. These results should be compared 
with numerical simulations [HHI- An exception to the above argument is the case of 
the IR attractor, namely the (3 function has a pole. In this case, we have not yet 
found the consistent solution. 

In section 7, we analyze the modified MA gauge. The MA gauge looks quite 
different reflecting a fact that the propagators are gauge dependent quantities. We 
show that the off-diagonal gluon and off-diagonal ghost propagators converge to Q^- 
independent constants in the IR limit. This implies that the off-diagonal modes 
become massive as is expected from the Abelian dominance jSUHZOJI^. The Abelian 
dominance was predicted in [70] and subsequently confirmed by numerical simulations 
on a lattice [TJ. On the other hand, independent analytical investigations suggest 
that the vacuum condensates of mass dimension 2 may be an origin of dynamical mass 
generation for off-diagonal modes, see [^n|l4ip42|E3] for the BRST non-invariant ghost 
condensate and |HH|IH7j for the BRST-invariant gluon-ghost condensate. However, the 
precise relationship between the condensate and the IR asymptotic behavior is not 
yet fully understood. A remarkable result in this article is that the diagonal gluon 
propagator can be of massive type and that only the diagonal ghost propagator can 
be short-range, since only the diagonal ghost does not obey superconvergence. The 
implication of this result to color confinement will be discussed in this article. Our 
results can be compared with the recent numerical results [7^ l73j . 
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The final section is devoted to conclusion and discussion. Some of the details of 
the calculations are given in Appendix A. 

2 The generalized Lorentz gauge and Maximal Abelian 
gauge 

In this section we introduce the generalized Lorentz gauge and the modified Maximal 
Abelian (MA) gauge which are used in this article for fixing the gauge of the SU{Nc) 
Yang-Mills theory. 

2.1 Generalized Lorentz gauge 

First, we modify the Lorentz gauge into more general form so that the usual Lorentz 
gauge is obtained as a special case. The gauge fixing (GF) term and the associated 
Faddeev-Popov (FP) ghost term for the generalized Lorentz gauge (so-called the 
Curci-Ferrari gauge jSSj) is combined into a compact form with two gauge- fixing 
parameters a and a' (See |2nilZlllZ3l2ZllZnil2Hll2nilZZllZHllZnilHniinnilHII): 

^GF+FP = ^^B^B (^1^ ■ - f ■ ^) + l^i^ ■ ^, (2.1) 

where ^ and ^ are the non- Abelian gauge field, the Nakanishi-Lautrup (NL) 

auxiliary field, the Faddeev-Popov (FP) ghost and antighost fields respectively, and 
(5b and Sb are the Becchi-Rouet-Stora-Tyutin (BRST) jH21 and anti-BRST transfor- 
mations [HB]- We use the notation 

F-G:=F^G^, F^:=F-F, {F x G)^ := f^'^^ F^G^ , (2.2) 

with /^^c* \yQYn.g the structure constants of the Lie algebra 5f of the gauge group G 
and A, 5, C = 1, ■ ■ ■ , A*"^ — 1 for G = SU (Nc). We adopt the convention in which the 
ghost and antighost fields are Hermit ian: 

^^(x) = ^(x), ^^(x) = ^(x), (2.3) 

and they are independent from each other. The pure Yang-Mills Lagrangian and the 
GF+FP term are separately invariant under the BRST and anti-BRST transforma- 
tions. 

By performing the BRST and anti-BRST transformations explicitly, we obtain 

^GP+pp = ^^^^ ■ ^ + ^ ■ df,£/^ - |^(^ X i^) ■ 

+ i^i ■ d^,^^\s^Y€ + ^/(^ X '^) ■ (^ X ^) (2.4) 
8 

= ^^^^ ■ ^ + ^ ■ d^£^^ - ^gi't^ xi^)-^ 

+ i^i ■ d^^^{£^Y€ + ^/(<^ X i^i) ■ (^ X z^). (2.5) 
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The GF+FP term includes the ghost self-interaction whose strength is proportional 
to the parameter a. In the limit a — > 0, the above GF+FP term reduces to the usual 
form of the Lorentz gauge which is found in the most textbooks of quantum field 
theory, 

^GF+FP =^<5b5b ■ + ■ ^ (2.6a) 

= - 2(5b ■ d^s^^ + —SS ■ (2.6b) 

■ d^^^{s^\€ + m ■ d^s^^ + ■ ^. (2.6c) 



The Landau gauge corresponds to a' = in the usual Lorentz gauge. 

For later convenience, we introduce another parameterization of two gauge pa- 
rameters: 

A:=« + a', e:=^(=^), (2-7) 
a + a' V IX) 

to rewrite the GF+FP term into 

A 
2 



+ ■ d^^^\£^\€ + X ^) ■ (<^ X ^) (2.8) 

+ Vi ■ d^^^\s^\€ + Y^'C^ X z^) ■ (^ X z<^). (2.9) 

The GF+FP term ()2.9p has the reflection symmetry: the Lagrangian is invariant 
under the reflection with respect to = 1/2 and the simultaneous exchange between 
the ghost and antighost fields, i.e., ^^\ — Therefore, the Yang-Mills 

theory with = GF+FP term is equivalent to that with = 1 GF+FP term. The 
= case is usually used in the textbooks of quantum field theory. At ^ = 1/2, 
the Lagrangian is invariant under the exchange of ghost and antighost. For ^ = 1/2, 
moreover, the GF+FP term has a hidden super symmetry 05^(412), see |2H1- The 
Landau gauge is reached by various limits, e.g., A ^ for ,^ = 1/2 (a' = 0, a — > 0), or 
A ^ for ,^ = (a = 0, a' — > 0) in the usual Lorentz gauge. Therefore, the Landau 
gauge can be regarded as a common point for the three invariant subspaces specified 
by ^ = 0, 1/2, 1, see Fig. 4 of [021 Appendix A for more details. 

2.2 Modified Maximal Abelian gauge 

The GF+FP term in the modified Maximal Abelian (MA) gauge |22] is given by 



S'gf+fp = J d'^xid^dB 



(2.10) 



where the respective field = .^^^ is decomposed into the diagonal 

and off-diagonal components (A = 1, ■ ■ ■ , — 1 for SU{Nc) gauge group), 

$(x) = $^(x)T^ = ^\x)T + <l>"(x)T", (2.11) 
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and the index i runs over the A'c — 1 diagonal components corresponding to the 
maximal torus group [7(1)^'="^ of the gauge group SU{Nc) and the index a runs 
over the remaining A^^ — off-diagonal components for SU{Nc)/U{l)'^''~^ . The 
conventional choices are a = l,2;i = 3 for SU(2) and a = 1, 2, 4, 5, 6, 7; z = 3,8 for 
SU{3). 

By performing the BRST transformation explicitly, we obtain 



+ -g'^pb^fd^C''&C'C^ + -g-^f^b'^f'^^d^b^c^t^d ^ ^gf^bc-^b^ja^c 

8 4 2 

-agr^HB^C^C' + '^g^ pb^ ^cdi^a^b^c^dy (2.12) 

In particular, the G = SU (2) case is greatly simplified as 

S'gf+fp = J d^xl^B''D^[a]'''A^^' + ^B'^B'^ 

+iC''D^[a]^'D^'[aY''C'' - ig^e^'^e^^C^C'' A^"" A^ 
+iC''ge''\D^[af''Al)C'' 

-age'^^iB'^C^C^ + ^g'^e^'^e^'^C'C^CC'^'^, (2.13) 
which reduces after integrating out the NL field i?" to 

S'gf+fp = j d'x[- ^{D^ia^ A^'f 

+iC''D^[a]'"'D^'[aY^C'' - ig'^e^'^e'^^CC^ A^^^A^ 

_^^g2^ab^cd^a^b(jc^dy (2.14) 

In order to completely fix the gauge, we need the GF+FP term also for the 
diagonal field. A simple choice is the Lorentz gauge for the diagonal part: 



Sqf%fp — J d X idji 



^^|9m;(x) + ^5^(x) 



=S*a^a; + ^{By + iC'd^C' + iC'd^'{gf^''AlC'). (2.15) 

There are many other choices for the MA gauge. A characteristic feature of the 
MA gauge is the explicit breaking of global color symmetry in addition to the local 
gauge symmetry and only the global U{1) symmetry remains unbroken after fixing 
the residual local U{1) symmetry. For more details of the symmetry in the most 
general MA gauge, see [HTj . 
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3 Analyticity, multiplicative renormalizability and 
RG equation 

In this section we summarize for reader's convenience the basic ingredients which are 
necessary in the following argument, since they are scattered in various references. 

3.1 Renormalization transformation and asymptotic behav- 
ior 

We define the renormalization transformation of the field $ by multiplying the field 
by a finite positive factor Z 

«l>(x,^7'^/x'^A') = 4^'<^>(x,/,/i^A) (s.i) 

where g is the coupling constant, /i the renormalization scale, and A is the gauge fixing 
parameter. Such a transformation is called the equivalent transformation. Two sets 
of parameters [g^ /i. A) and [g' , /i'. A') are called equivalent parameter set and are 
denoted as {g,fi,\) ~ {g',fi',\') if the corresponding field operators are related by 
the renormalization transformation ()3.1|) . 

We define the effective coupling or invariant charge 

Q = Q(A;^^?^/i^) = Q(^,/) (3.2) 

as a dimensionless invariant of the RG which satisfies 

Q{k'^,g^,f?) = g^ at k"^ = f?. (3.3) 

Here the invariance under the RG means that the relation 

Q(A;^^?'^/i'^) = Q(A;^/,/i^) (3.4) 

holds for equivalent parameter set {g,fi) ~ {g',fi'). The effective coupling is usually 
defined in terms of time-ordered functions and an explicit example will be given later. 
We have required that the effective coupling or invariant charge should not depend 
on the gauge parameter A. 

By putting A;^ = //'^ in eq. ()3.4p . we obtain 

Q{^^'^g'^^l") = Q{^^'^g^^^') = g'^ (3.5) 

where we have used eq. ()3.3|1 in the last equality. Then the definition ()3.2j) implies 

«'^^sg,.^). (3.6) 

The propagator D of the field $ is transformed like 

D{e,g",fi'\X') = Z^D{k\g\fi',\). (3.7) 
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Moreover, we introduce the dimensionless function, say, the dressed function (or form 
factor) R by 

i?f^,/,A) ■.= -eD{e,g^fi',\). (3.8) 



9 ' 



Then R is transformed hke 



R(^^,g'\X')=Z^R(^^,g',xY (3.9) 



By putting k = /i in this equation, we obtain 

,'2 



R{l,g",X')=Z^Ri^,g\X\, (3.10) 



and hence Z$ is written as 



Z^ = R-'(^^,g\X^R{l,g",X'). (3.11) 



Thus the transformation ()3.9p of R is rewritten into 



R { ^, a") i? ("1^, X'j = R (1, A') R r-, g\X]. (3.12) 



Finally, the gauge parameter transforms hke^ 

X' = X\g'\X) = XZx, (3.13a) 
and in particular, for A = 0, ,^ = or ,^ = 1 

A' = X\g", A) = XZ,^ = XR~' (^^, g\ X^ R (l, g", A') . (3.13b) 

The conventional normalization condition for the propagator: 

-k^D{k^,g^,fj.^,X) = 1 at /c^ = /^^ < (3.14) 
leads to a normalization: 

i?(l,/,A) = l. (3.15) 

Thus, under the conventional normalization condition, the renormalization transfor- 
mation of R is given by 



R ( ^, X]=R 9\ a) R A') , g- = Q g^ ] (3.16) 



Here the gauge parameter transforms in particular for A = 0, = or ,^ = 1 like 



A' = X'{g", A) = XZ,, = XR-' g\ A^ (3.17) 
with the normalization 

A'(/,A) = A. (3.18) 



^In the conventional Lorentz gauge (corresponding to ^ = 0), it is well known that Z\ =■ Zi;^ . 
However, this is not the case in the generalized Lorentz gauge ^ ^ 0, 1 except for the Landau 
gauge A = 0. In fact, the one-loop calculation |^ in the dimensional rcgularization leads to 



Zx = Z^ + Ae(l - """f^ where e := 2 - D/2. 
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3.2 Non-perturbative multiplicative renormalizability and RG 
equation 

We require the non-perturbative multiplicative renormalizability. By introducing the 
renormahzation factor Z$ for the field $ with the renormahzation scale fi and the 
ultraviolet cutoff A, the bare (unrenormalized) n-point function f g"^ is related to the 
renormalized n-point function T^^'> as 

fl;'\p^,■■■ ,p„.;A^aoAo) = Z^ifx^A'')-^^'f^^\p,,--- ^Pn^lJ^^a^X), (3.19) 

where is the bare coupling and a := g^/4:7r is the renormalized coupling constant. 

The differential form of the renormalization transformation for the Green functions 
is the renormalization group (RG) equation. The trivial relation 



^rJ"^(Pi, ■ • • ^Pn, A^ao, Ao) = 



(3.20) 



implies that the renormalized n-point function F^") obeys the RG equation:'' 

f(")(pi,-- - ,p„;/i^«,A) = 0, (3.21) 



where we have defined the RG functions by® 

da{n) 1 91nZ$(/i) 1 d\nZx{fi) 



(3.25) 



It is shown jHl] that the general solution of the RG equation ()3.21|1 is given by 
f(")(spi,-- - ,sp„;/i^a,A) 

-n dt'-f<s,{a{t')) 



/ ""exp 



-n / da 



f^^\p^,--- ,Pr,.-^\a{t)Xt)) 
f(")(pi,--- ,p„;/i2,a(t),A(t)), 



t :-- 



m d\ 



In s, 



(3.26) 
(3.27) 



^For the generahzed Lorentz gauge, we assume ^ is fixed. In fact, this is possible for ^ = 0, 1/2, 1, 
since the relevant Lagrangian has distinct global symmetry and the RG flows are restricted into the 
subspace specified by^ = 0,l/2,l respectively, see 

^The renormalization relation gu = Z~^g implies the mutual relationship: 



I3{g) 1 dg 1 dg^ p{a) 



Therefore, the (3 function can be calculated by differentiating Zg by fi. 



9 In 9 



dh-iZa 



f3{a) = -2afi 



dlnZg 
9/1 



Here we have used the relation X = XqZ^ to yield 

d\ , dlnX , dlnXoZZ^ 



-A/i — = -2A7A. 

a/i 



(3.22) 



(3.23) 



(3.24) 
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where we have taken into account the fact that dim[f ("^] = 4 — n for the field $ with 
dim[$] = 1, and a{t) and X{t) are the running couphng constant and the running 
gauge parameter. 

In particular, the unrenormalized (transverse) gluon and ghost form factors Fq, Gq 
are related to the renormalized ones Fr, Gr as 



Fo(/, A^ ao, Ao) =^3(/i^ ^'^)Fr{p^, /i^ a, A), 
Go(p^ A^ ao, Ao) =^3(/i^ ^^)Gr{p^, /i^ a. A). 



The renormalized form factor Fr obey the RG equation: 

FR(p\^^\g\\)=Q. 
The general solution is given by 



(3.28) 
(3.29) 



(3.30) 



Fr{sp ,a. A) = exp 



F«(p^/i^a(t),A(t)). 



(3.31) 



In addition, for the general form of the full gluon propagator given by 



F(y 

p2 



p2 



+ A 



F^{p^)PtiPv 

9 9 ' 



(3.32) 



the unrenormalized longitudinal gluon form factor Fq is related to the renormalized 



one F^ as 



F^{p\ K\ ao, Ao) =Z^\f^', A')Z,ifi^ A^)F^(p^ /x^ «, A). (3.33) 
The renormalized longitudinal form factor F^ obeys the RG equation: 



9 ^/ 9x 



2A7a(/)^ + 27.^(/)-27a(/) 



F^{p\fi',g',\)=0. (3.34) 



The general solution is given by 



^_R(■5p^/i^«, A) = exp 



— W) — 



F^(p2;/i2,a(t),A(t)). (3.35) 



In the conventional Lorentz gauge, 7^/(0) = 7a(«) and the longitudinal form factor 
of the gluon does not run, i.e., Fj^ijP' , jj?, a. A) is a ^^-independent constant, as is well 
known. This is not the case in the generalized Lorentz gauge. 

It is well known that the Yang-Mills theory is multiplicatively renormalizable at 
least in the perturbation theory. The RG functions can be reliably calculated by 
perturbation theory in the UV region by virtue of asymptotic freedom 



/^(") = -^"^ ' 7<i'(") 

ZTT 



An 



(3.36) 



13 



Therefore, the Yang-Mills theory has an UV fixed point at a = 0. The UV asymptotic 
form of the n-point function is given by 



:''0 



= ''~"(f[iy) ''°^("HPl,■■■,Pn;/i^a(t),A(^)). (3.37) 
Thus the UV behavior of the propagator (p^ — > oo) is given by 



1 f a{p) 

VP) ~ 

for the running coupling constant 



^0 



3.3 Renormalization group on the complex plane 

The equation ()3.12j) or ()3.2H) is valid also for complex k^. Then the functions D and 
R satisfy the same RG equation: 



r(^^, 9^X^=0, (3.40) 



with the same coefficients, i.e., RG functions as those of ()3.2H) at every point k"^ of 
the cut k'^ plane (i.e., complex k^ plane except for the positive real axis). 

We consider the propagator along the ray in the cut k"^ plane (See Figure |21): 

= ^\k^\e'^ = \k^\e'^^+y^\ -n < ^ < n, (3.41) 

where the angle ip is measured from the negative real axis A;^ < in the direction 
indicated by the arrow in Figure |21 Setting /i'^ = — the relation (j3.12|l yields 

R (^^, x^=R (^^, a) r a) R-' (l, a) , f = Q (^^, / 

(3.42) 

which is rewritten as 

R (^^, X^=R (^^, a) R (e^^ g. A) R~' (l, ^^ A) , f = Q=Q ^7^) . 

(3.43) 

In addition, the gauge parameter is transformed in particular for A = 0, = or 
^ = 1 as 



A = XR~' a] R a) . (3.44) 
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9^"'---. 



Imk 



X X X X X X X >- 2 

(kSo) ^ek 

Minkowski 



(k^<0) 

Euclid 



Figure 2: The ray with angle ip measured from the negative real axis /c^ < on the 
complex k"^ plane with singularities on the positive real axis. 



Using the asymptotic behavior of the structure function for g'^ +0 
R{e, g\ fi^ X) = R f ^, \] = 1 + j^{\)g' In ^ + 0{g'), 



and the continuation, 



we obtain 



(e*^ Q, a) - 1 + 7o*(A)^V + 



\k2\ 



lim R fe^^, Q, A) = lim R (e'^, Q, X) = 1. 



(3.45) 



(3.46) 



(3.47) 



Now, it is important to recall the asymptotic formula in Euclidean region fc^ < 
which is obtained from the general solution of the RG equation ()3.40j) : 



c^(A) = {g'mr^ow/^o 



2/ c^x^ 



exp 







C$(A) In 



> 0. 



-7,?(A)//3o 



(3.48) 
(3.49) 



In the conventional normalization i?(l,Q,A) = 1. In other cases, the leading be- 
havior of the coefficient i?(l, Q, A) for g"^ ^ +0 or m := k'^/fi^ — > oo is determined 
by the asymptotic form of R{l,g^,X) for small g. The behavior of the exponen- 
tial factor follows from the behavior of Q and the ratio 7$//5. Here the integrand 
7$(a;)//3(x) is not integrable at x = in QCD and therefore the singular part is sep- 
arated 7$(x)//9(x) = ■Jq / PoX~^ + r$(x) with r<i,(x) being the remainder integrable at 
X = 0. The coefficient C$ is positive. See Appendix A. 
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Thus we find the asymptotic form for |A;^| ^ oo along the ray ()3.4H) : 

■ \k' 



.2|\ -7f(A)//3o 



R{k\ g\ ii\ A) =Ras{k\ g\ //^ A) = C$(A) In 



D{k\g\^\\) ^D,s{k\g\^\\) = -C^{X)k-' In 



|^2|\ -7o*(A)//3o 



3.4 Analyticity 



dxT{x) 



> 0. 




(3.50) 

(3.51) 
(3.52) 



Figure 3: An integration contour C on the complex k"^ plane with singularities on the 
positive real axis. 



If the analytic function f{z) vanishes along any ray in the cut complex plane as 
\z\ — s> oo, then it has the dispersion relation as follows. 

Suppose that the complex function f{z) is analytic in the whole complex z := k^ 
plane except for the positive real axis z > 0. Then we choose a closed loop C in the 
complex fc^ plane such that the complex function f{z) is analytic inside and on the 
closed loop C as in Figure El For a reference point k"^ inside C, the Cauchy integral 
formula tells us that for the closed contour C of integration 



2Txi Jc 



dz 



A;2 



1 1^ , J\z + ie)- f{z-ie) ^ 1 / 
liri Jo z — k^ 2Tii J\z\=R 



(3.53) 

where we have separated the integral into two pieces; one is the contribution from 
the paths above and below the positive real axis; the other is the contribution to the 
integral over the circle of radius R. 
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First, making use of the fact that /(/c^) is real for k"^ < s„nn (at least for the space- 
like region k'^ < 0), we can use the Schwartz reflection principle f{z — ie) = f*{z + ie) 
for Rez > Smin and hence f{z + ie) — f{z — ie) = f{z + ie) — f*{z + ie) = 2ilmf{z + ie). 
Thus the first term in RHS of ()3.54j) picks up the discontinuity (the imaginary part) of 
the f{z) function along the positive real axis. Next, the circle \z\ = Ris parameterized 
&s z = Re^^: 

ne) = i + f jg^f/(^f:). (3.55) 



TT JO z — k'^ Jo 2tt i?e*^ — /c^ 

The second term in RHS of this equation has the upper bound 

de Re'^fiRe'^) f^- dO R\f(Re'^)\ ,..„,flM R 

Now we wish to put the radius R of the circle to infinity. Thus, if \f{z)\ ^ as 
\z\ — * oo, then the second term has no contribution as R ^ oo and we have the 
dispersion relation:^ 



This is also written as 



/(P) = r ds^^^^, p{s) := -Imf{s + ^0), (3.58) 

where k"^ is an arbitrary point in the complex plane except for the positive real axis. 
Note that f{z) need not to go to zero in the power-like, i.e., f{z) ~ l/|z|^(£: > 0) as 
is usually assumed. In fact, a weaker decay f{z) ~ (log |z|)~^(£: > 0) is also allowed. 
Incidentally, since the Cauchy integral formula is generalized as 

the similar argument as above leads to 

/""(*^)-'/;^»(j^- 

Note that the dispersion relations ()3.57|) or ()3.58j) is derived from 

i) the analyticity of f{z) except for the positive real axis, 

ii) the asymptotic vanishing of f{z) as \z\ oo, 

iii) the real-valuedness of f{z) on the real axis z < Smin (at least for the negative 
real axis). 

In the next section, we equate the dispersion relation with the spectral representation 
to extract the information on the spectrum using the asymptotic behavior of the 
propagator and the dressed functions. 

^In general, the boundary value of the analytic function is defined as the Sato hyperfunction, 
which is an extension of the Schwartz distribution. In applying this relation to the propagator and 
the form factor later, we require the positivity in Euclidean region, i.e. on the negative real axis of 
the complex fc^ plane. Then they become a positive measure and can not be a distribution. 
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3.5 Minkowski region 

In general, the propagators are tempered distributions in the Minkowski region > 0. 
However, if distribution-type singularities are present only below a finite momentum 
value K, the propagators may be treated as ordinary functions for k"^ > K"^ [IH]- In 
this case, therefore, the asymptotic formulas derived above are also applied to the 
large Minkowski momenta with ip = ivr. 

If the distribution-type singularities occur for arbitrarily large momenta, the 
asymptotic behavior of the propagators should be studied in terms of the average 
values which are obtained by smearing with suitable test functions. It is shown 
that the averaged gluon propagator in the Minkowski region has again the same 
asymptotic form as given above, independent of the chosen test function. 

The asymptotic behavior for the weight function of the spectral representation is 
derived as follows. The weight function has the expression: 

7rp{k^) = lmD{P + i€) = -k-^lmR{k^ + ie), P > 0, (3.61) 

where we have used the result of the manipulation: lm[R{k'^ + ie)] = —lm[{k'^ + 
ie)D{k^ + ie)] = -kHmD{e + ie) - eReD{k^ + ie). Setting //'^ = in (jTT^ . we 
have 

R x]=R a) R {-l,f, A) (l, f, A) , f = Q f^, gA , 



(3.62) 



which leads to the imaginary part for fi < 



7rp{e) = -k-HmR{e + ie) = -k-^R g\ A j ImR [-l,g\ A) R'^ (l, ^^ A) , 

(3.63) 

where R (j^, g'^, A^ and R (^1, (7^, A^ are real for k'^ > 0. From ()3.46|) . we obtain 

Imi?(-1,Q,A) ='iT^UX)Q+Q^h, (3.64) 

with the remainder h vanishing for Q — > 0. It is shown by one-loop resumed pertur- 
bation theory that 

Q-47r|/3o-^| (^In^^ . (3.65) 

Thus we obtain the asymptotic form 

p(/^;^5(^/i^ A) =Pa3(/^;^5(^/i^ A) as /c^ 00, (3.66) 



pUe, g\ /x^ A) = - (A:2)-So*(A) QRas [j^y g\ A j (3.67) 
= -47r(A;^)-So*(A)|/3o"'|C*(A) (^In^j , (3.68) 



18 



where it is important to remember that the anomalous dimension 7$ is gauge depen- 
dent and so is the asymptotic form. 

If the weight function has distribution-type singularities at arbitrarily large mo- 
menta, we consider the asymptotic behavior of the average values. Consequently, the 
same asymptotic behavior is obtained j46j . 

It should be remarked that the sign of the asymptotic discontinuity pas is deter- 
mined by the ratio jo/Po, i-^-, the sign of p is the opposite of'jo/Po. It is negative for 
10/ > 0. 



4 Naive derivation of superconvergence relation 

Now we can derive the superconvergence relations in the generalized Lorentz gauge. 
Our results reproduce the well-known superconvergence relation in the conventional 
Lorentz gauge. The details of calculations are given in Appendix A. 



4.1 Gluon 

For 'Jq/Po > 0, the renormalized gluon structure function Dt has the property that 
the complex function defined by 



Dr{k\p\gl\a) + ^ (4.1) 

vanishes faster than indeed it behaves like {k'^)~^CA (in % /^o large 

\k'^\ in all directions of the complex k"^ plane. For ■Jq/Po > 0, moreover, the renormal- 
ized dressed function R = —k^Dx has the property that the complex function defined 
by 

R[k\p\glAR)-Y^ (4.2) 

vanishes like Ca — >■ as \k'^\ — > cxD, in all directions of the complex k'^ 
plane. By applying the formula (j3.57|) to the function ()4.1|) . therefore, we can write 
an unsubtracted dispersion relation 

DAk ,/i ,gR,XR) = / dp , (4.3) 

^0 P rv 

where we have used a mathematical identity: 

1 



P—-i'k5{p'). (4.4) 



p"^ + ie p2 

In the similar way, the formula (j3.57|) applied to the function ()4.2p yields 
2 2 2 ^ ^ ^ /"^ , 2PV(P^/^^^R,Ai^) 



A* JO — k-^ 
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Therefore there are two alternative forms for the spectral representation for D: 

n (h^ n2 \ \ r ^^2 piP^f^^9RAR) A/A. 1 ^ p{p^ , fi^ , Q^, Xr) 

Dr{k,^^,g^,Xn)-l dp ^-^^ --^+^70 ^^^F 

(4.6) 



This is possible if and only if the relation 

A 



;>oo 

/ c^pV(p^ 

Jo 



A* 



(4.7) 



is satisfied for arbitrary fc^. For / (3q > 0, thus, we have the superconvergence 
relation in the generalized Lorentz gauge with an initial gauge- fixing parameter A: 

POO \ 

/ dep{e,fi',gR,XR) = -. (4.8) 
Jo A* 

This relation reduces to the well-known superconvergence relation in the Landau 
gauge 

;>oo 

/ dep{e,p^gR,0)=0. (4.9) 
Jo 

Here A* depends on ^ in the generalized Lorentz gauge: A,, = 13/3 in the conventional 
Lorentz gauge ^ = and A* = ^ ^ ■'-Z^' ~ T| ^ ^ see section V 

of [nn] and Appendix A. 

Irrespective of the gauge parameter, therefore, the gluon form factor vanishes in 
the IR limit due to superconvergence relation ()4.8|1 : 



A 

R{0, p\ gl Xr) = -- dp^p{p\ p\ gl, Xr) = 0, (4.10) 
A* J-o 

and hence the gluon propagator and the form factor has the gauge-parameter inde- 
pendent IR behavior for the small momenta: 

F(Q2, p', gl Xr) = 0(g2), D{Q\ /i^, g^ Xr) = const. + O(Q^). (4.11) 

These behavior are independent from the gauge parameter. 



4.2 Ghost 

However, the situation is quite different for the ghost. We can repeat the same steps 
as in the gluon case. Then we obtain the following results which agree with those 
of jSni obtained in the conventional Landau gauge. In the arbitrary gauge, the ghost 
has an unsubtracted dispersion relation for the renormalized ghost propagator App: 

AMk\ ,\ 9l 0) = r (4.12) 

JO p — /e 

However, the form factor Gpp := —k'^App has an unsubtracted dispersion relation 
only in the Landau gauge A = 0. 

Gpp{k ,p ,gR,0) =- J^^ dp ^— p . (4.13) 
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Hence, the super convergence relation for the ghost holds only in the Landau gauge 
A = 

0= depFpik^fi^glO). (4.14) 
Jo 

For A 7^ 0, however, the superconvergence relation does not hold, since the unsub- 
tracted dispersion relation exists only for the propagator, not for the form factor: 

Appik , /i ,gji, Xr) = dp ^2rrp = ~~gi'^PPy^ ^1^ i9r, ^r)- (4.15) 

If this result is true, then the superconvergence for the ghost holds only in the Landau 
gauge, in sharp contrast with the gluon case. There is a subtle point which may 
invalidate the above derivation. We reconsider the ghost in the next section. 

5 Superconvergence relation and spectrum 

In this section, we combine the superconvergence relation with the spectral represen- 
tation to discuss the spectrum of the Yang-Mills theory. 

5.1 Spectral representation 

From the assumptions: 1) Poincare group structure (representation), 2) spectral con- 
dition, 3) completeness condition, the spectral representation (Umezawa-Kamefuchi- 
Kallen-Lehmann representation) |HS1 for the gluon two-point function follows: 

= f (5.1) 

where p((T^) is the spectral function. Due to a mathematical identity 

= P- + i7i6(x), (5.2) 
X — le X 

the spectral representation 1)5.11) of the propagator Dxik"^) = ReDxik'^) +iImDT{k'^) 
is decomposed into two parts: 

ImDrik^) = 7rp{k^), (5.3) 

and 

ReDrik') = / da'pia')P^—. (5.4) 

Jo (7 K 

Thus we obtain the same relation as the dispersion relation: 

ReDTik^) = P da^pia^)^—, pip^) = -ImDrip^). (5.5) 

However, it should be remarked that we have not used the analyticity for deriving 
this representation. 
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If the state space has positive definite inner product, the spectral function must 
be non-negative: 

p(p^) > 0. (5.6) 

In the state space with the indefinite inner product, the complete set {\n)} can not 
be spanned by the eigenstates of alone in general and the metric 

{n\n') = r]nn' (5.7) 

can not be diagonalized, and the completeness condition must be replaced by 

^ = T.\^)Vnn'{n'\. (5.8) 

n,n' 

In Yang-Mills theory, fortunately, all |n) can be taken to be eigenstates of and 
the metric can be diagonalized into the form = ±1. In this case, the positivity of 
the spectral function does not hold, since both positive and negative terms contribute 
to the sum defining the spectral function. This fact plays the important role in the 
superconvergence relation in Yang-Mills theory. 



5.2 IR behavior of the renormahzed gluon propagator in Eu- 
cUdean region 

The spectral representation of the renormalized gluon structure function reads 

Drik ,/i ,^r,Ar) = / dp , (5.9) 

^0 P rv 

P(P^ 9% Ar) =n~'^ImDT{k^ + ie, /i^ g\, Xr). (5.10) 

The renormalized version of p is a function of p^, /i, qr, Xr where fi is the renormaliza- 
tion point. Note that this representation holds for arbitrary point k"^ in the complex 
k"^ plane except for the positive real axis k^ > 0. It was shown that the gluon structure 
function has the superconvergence relation: 

POO \ 

/ dp^p{p^,i?,g\,XR) = (5.11) 

in the parameter region satisfying 'Jq/Po > 0. In this paper we consider only the case 
/5o < in pure Yang-Mills theory corresponding to asymptotic freedom. Defining the 
dimensionless dressed function R by 

R{e) := -eOrik^), (5.12) 

we can write the dispersion and spectral representation of the dressed function: 

R{k ,1^ ,9r,Xr) = Y - dp ^2 _ fc2 • (^-^3) 

Putting the reference point /c^ on the Euclidean region, i.e., = —k"^ > 0, the 
gluon structure function has the representation: 

Dt{-Q ,/x ,9r,^r) = dp ^ . (5.14) 
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In fact, Dt is not singular in the Euclidean region > 0. Defining the dimensionless 
gluon form factor F by 

F{Q') := R{-Q') = Q^Dt{-Q'), (5.15) 

we can write the dispersion and spectral representation of the form factor in the 
Euclidean region Q"^ = —k^ > 0: 

F{Q ,1^ ,9r,Xr) = ^ - dp . (5.16) 



A, Jo " p2 + g2 

Confinement is considered to be a strong coupling problem and to be related to the 
IR region. Therefore, we focus on the IR behavior of the gluon form factor. We find 

that the gluon form factor vanishes in the IR limit of the Euclidean region +0: 

A 

F(g2 = 0, /i^ gl Xn) = -- dp'p^p^ g^ Xr) = 0, (5.17) 

A* JO 

due to the superconvergence relation (j5.11|l . Moreover, the spectral representation 
has the formal power series for small Q^, 



roc 

F{Q\n\glAR)=T.(-ir\Q'T dp' 

1 JO 
n=l 



P{p'^,f^^,9R,^R) 



(p2)n 

p2 ^^'^'jo (p2)2 



__Q. I ^^.pjAjA9l^ _ ^^, p{p\,\glAR) ^ ^^^e^ 

(5.18) 



while the gluon structure function has the IR behavior 
Dr{-Q\p\ gl Xr) = £ (-l)'^(g^)" f dp'^^^jf^^^ 



n=0 



(p2)n 

dp^PjllllAM _ g2 r ^^, p{p\,\9i,XR) ^ 

Jo no"' r 



p2 Jo ip^ 

(5.19) 

In the Euclidean IR limit i 0, the gluon structure function converges to a constant: 



n 2 2 X ^ r ^ 2 p(p^/^^^l^Afl) 

^r(0,/i ,gR,XR) = / dp , (5.20) 

JO p^ 



provided that all the integrals are convergent. If this constant is non- vanishing and 
finite, we have k = 1/2. 

It is expected that each integral appearing in the coefficient of the power series 
()5.18|) and ()5.19|) is positive and the series is alternating. In particular, the constant 
Dt{0) ()5.20|) is considered to be positive. The integral is convergent for large p' due 
to the asymptotic form ()3.68|) . For large p^, the integrand p/p"^ is small and negative, 
while it is large positive for small and it does not diverge due to the existence of 
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mass gap. Therefore, the total integration is expected to be positive. This conclusion 
is a gauge independent result, i.e., independent of the choice of gauge parameter. 
If the constant ()5.20|) happens to vanish, 

r-^^, p{p\,\ghXn) (5.21) 

JO 

then the k = 1 case occurs, 

F{Q', 9l, Xn) = ~{Q'? f dp^P^^l^^^l^ + OiiQr), (5.22) 

and 

Dr{-Q^f.\9lXn) = -Q' f dp'^^^^^^^^^ + 0{{Qr). (5.23) 

At first sight, this case is similar to the Gribov case. However, there is a crucial 
difference and this case is unlikely to be realized as argued in the next subsection. 

The above arguments hold as far as the IR divergence does not occur. In other 
words, the integrals for the spectral function are assumed to be finite. This issue is 
discussed below. 



5.3 Moments of spectral function and mass gap 

Now we consider the implication of the above result to the spectrum. To this end, 
we consider the moments of the spectral function: 

rR 

J dss^p{s) (£ = 0,±1,±2,--.). (5.24) 

This integral is convergent for large i?, if ^ < 1 due to the asymptotic form ()H.fi8j) of 
p{s) for large s = k"^. Here any negative value of i is allowed. In what follows, we 
pay attention to the convergence for the limit e — >■ 0. 

To obtain the moment, we apply the Cauchy theorem to the function f{z) which is 
analytic inside and on the contour C in the complex plane shown in Fig. |3] avoiding 
the positive real axis and the origin. Then the function f{z) satisfies 

i dz z^f{z) = 0. (5.25) 

Here we separate this integral into three pieces; one is the contribution from the 
paths above and below the positive real axis [e, R] which pick up the discontinuity 
(the imaginary part) of the function f{z) along the positive real axis; the others are 

"'^"Now we consider the implication of the above result to the spectrum. If the gluon spectral 
function has a massless pole and the single particle contribution is explicitly separated as plp^) — 
Zd{p^) + • • • , then Dt{0) given by H5.20() diverges. Therefore, finiteness of the renormalizcd gluon 
propagator at = excludes the massless gluon pole, suggesting the existence of mass gap. In 
this case, all the coefficients of the power series in are finite due to the asymptotic form H3.68|l 
of p for large p^. Consequently, the IR power series expansions (|5.18() and H5.19|l arc well-defined. 
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the contribution to the integral over the circle of radius 1 and the contribution 
to the integral over the circle of radius e ^ 1. Hence ()5.25j) is decomposed as 



ds s^-Imfis + zO) =— / dz z^fiz) - — i dz z^f(z) (5.26) 
Je n 2m J\z\=e 2m J\z\=r 

/•2vr AO /•27r Afi 

Jo 2tt Jo 2tx 



(5.27) 

For I = —1, we apply ()5.27|1 to the gluon structure function / = Dt'- 

pR 1 /•27r Af) 1-2-K Af) 

/ ds s-^-ImDris + %{))= — Driee^') - — DriRe''). (5.28) 
Je 71 Jo 2tx Jo 2tx 

In the region where the super convergence relation holds, Dx{z) vanishes in the limit 
\z\ = R ^ oo. In the limit R ^ oo and e — 0, therefore, ()5.28|) reduces to 

/ dss-^p{s)= ds s-^-ImDris + iO) = — Dt{0) = Dt{0), (5.29) 
JO Jo vr JO 271 

where lim^^Q Dxiee^^) has a common value in the limit of approaching the origin 
= e — along the ray k"^ = |A;^|e*^(0 < 6 < 27t) except for the positive real axis 
9 = 0. This result ()5.29p is consistent with the previous result ()5.20|) . 

The finite Dt{0) is consistent with the spectral function with a massive pole (at 
p2 = m2): 

p{s) = Z6{s - M^) + p{s), (5.30) 
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where p{s) is the contribution of the continuous spectrum from more than two particle 
states, namely, p{s) has no support for s < (existence of the lower gap). This is 
because the finiteness of the IR limit -Dt(O) denoted by the expression 

DriO) = / ds s-'p{s) = YH+ ^"'p(^) (5-31) 

JO iW^ JO 

is interpreted as indicating the existence of a massive pole (M 7^ 0), provided that 
Z 0. In this sense, the finiteness of Dt{0) excludes the massless gluon pole: 

p{s) = Z6{s) + p{s). (5.32) 

The same argument holds for other moments of the spectral function appearing in 
the power-series expansions ()5.18|1 and ()5.19|) . Thus, we have shown that, for gluon 
with massive spectrum ()5.3()j) . then the power-series expansions ()5.18|1 and for 
the gluon propagator and the form factor could be well-defined for small Euclidean 
momenta Q^. 

Note that < -Dt(O) < 00 corresponds to k, = 1/2. This value for the IR 
exponent k seems to be supported by recent Monte Carlo simulations on a lattice ^\ . 
Quite recently, Bloch has proposed a truncation scheme respecting multiplicative 
renormalizability of the coupled Schwinger-Dyson equations for the gluon and ghost 
propagators by including all the diagrams [SB] to obtain the approximate solutions 
with the IR critical exponent k = 1/2. 

However, p{s) is not necessarily positive for s > 0. Therefore, a situation of 
vanishing Dt{0): ds s~^p{s) = -Dt(O) = happens to occur in principle. If 
this is the case, we must consider the next-to-leading term and then k = 1 will be 
realized (the apparent Gribov limit). However, we argue that Dt{0) = is unlikely 
to be realized based on the general properties of the spectral function. In fact, the 
structure function predicted by Gribov has the complex conjugate pair of poles 

This is incompatible with the dispersion relation (j5.14|l which requires the analyticity 
in the cut complex plane. 



5.4 IR behavior of the renormaUzed ghost propagator in Eu- 
chdean region 

For the Euclidean momenta Q, the ghost propagator has the spectral representation 

Afp(Q , /i , ^r, Ar) = dp p2j^Q2 = ql'^PPy^ ^1^ Afl)- (5.34) 

^^It is possible to consider the absence of one-particle pole Z = 0, i.e., p{s) = p{s) in Minkowski 
region. However, we suppose that the asymptotic fields exist, even if the corresponding quanta 
can not be observed (i.e., they are confined). Only in this formulation, we can adopt a sufficient 
condition for color confinement due to Kugo and Ojima. Therefore, we restrict our consideration to 
the case in which this framework works. 
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The propagator is not singular in Euclidean region > 0. Hence we consider the IR 
limit i 0. We see that the contribution from large to the integral is convergent 
due to the asymptotic form of the spectral function derived in Appendix A. For 
ghost without the super convergence relation, the IR limit of the ghost propagator is 
expected to be 



^Fp{Q i 0, /i , Aij) = / dp = oo. (5.35) 

JO p^ 



If (I5.35|) was finite, the ghost form factor had the similar form as the gluon case and 
the ghost form factor vanished in the IR limit: 

GMQ\ 9l Ah) = r dp^Pii^^liliAM + om - as - 0. 

JO p^ 

(5.36) 

Eq. ()5.35|1 suggests that the ghost propagator can not have the power-series expansion 
in around = in the Euclidean region. For the situation ()5.35p to be realized, 
the ghost can not have the spectral function of massive type. This result does not 
contradict with the consideration based on the RG equation given below and the 
approximate solutions of the truncated coupled SD equations. 
If the ghost has the massless spectrum of the usual type: 

PFp(p^/i^5'fl, A/j) = Zfp5{p^) + Pfp{p^), (5.37) 

with Zpp 7^ 0, then the ghost propagator reads 

Afp(<5 ,9r, Aij) = + dp p2^Q2 = qI'^fpKQ ,9r, Ar). 

(5.38) 

This again contradicts with the consideration based on the RG equation below. 

Thus we reconsider the derivation of the dispersion relation given above. If the 
spectral function p has singularities accumulating toward the origin p"^ = 0, we must 
replace an integration contour in Figure El by a new contour in Figured to avoid the 
origin. For a reference point k"^ inside a new contour C, the Cauchy integral formula 
yields 

" ^J(z + te) - f(z - ,e) 



2Tii Jo z — k'^ 

The similar argument leads to 

./,2n 1 fR , Imf(z + ie) f"^^ de Re'^fiRe'^) f^^ dO ee''^ f(ee'^) , 
^^^)=-/~ ~ P +70 2^ Re^^-k- -Jo 2^^^^- 



TT JO 
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The second term vanishes for the function f{z) which vanishes asymptotically as 
— oo and hence 

= i r iJj^l^ + r P4^. (5.42) 



7r JO 



z — k'^ Jo 27r k'^ — ee*^ 



The effect of the last term is classified into three cases. 

i) In the limit e — 0, the last term vanishes for the bounded function f{z) which 
is regular in the limit z ^ (massive case). Therefore, the above conclusion for the 
gluon with massive spectrum also holds for a new choice of the integration contour. 
This is an undesirable result. Therefore, we allow f{z) to be singular in the limit 
z ^ 0. For concreteness, we examine the asymptotic form f{z) ~ Cz^~^ for 1^1 -C 1. 
Nevertheless, vanishing of the last term also follows for the function which becomes 
singular in such a way that f{z) ~ Cz^~^ with positive real (3 > 0. This is also the 
case for a negative integer (3, since for sufficiently small e 

2^"fc^^^ ~ U^io 2^^ ""P 2mf3 ' ^' ' 

ii) For /3 = 0, i.e., f{z) ~ Cz~^ (the tree massless form), however, we have the 
additional term: 

iii) Finally, the additional term remains also for the singular function f{z) ~ 
Cz^~^ with negative and non-integer (3 (IR enhanced case). It is this case that 
corresponds to the ghost propagator in the Landau gauge, as argued below based on 
the RG equation. 



5.5 Renormalization condition 

This should be compared with the momentum subtraction scheme in perturbation 
theory where the renormalization condition fixes the gluon propagator to the tree-level 
one at a sufficiently large space-like renormalization point k"^ = fi^, i.e., /^^-(/i^, /i^, gj^, Xji) = 
This will lead to 

POO 

/ /^^ ^ oo, ^ 0, Xr A,) = 1. (5.45) 

Jo 

In the formal level, the spectral representation has the expansion for large 

^(Q^ 9% A/j) = / dp"^ pip^, /i^ gl, Xr) 

Jo + 

= E77«ir/ dp\p'rp{p^^^\glAR). (5.46) 
n=o I •'^ 

This power-series expansion in l/Q^ is not well-defined, since finiteness of the integral 
as the coefficient is not guaranteed for n > 1. 
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The relation ()5.45p seems to contradict with the super convergence relation. How- 
ever, the renormalization condition F{fi'^,fi'^) = 1 is too strong, since it specifies 
in advance a part of the non-perturbative dynamics which should be determined 
by solving the theory. Therefore, we do not require the renormalization condition 
= 1 even for large /i in this paper. Rather, we adopt the renormalization 
condition for any value of /i 

F{fi\fi^)G\fi\fi^) = l, (5.47) 

where G is the ghost form factor. A loophole of the above argument is to take into 
account the FP ghost in the super convergence relations. 
The super convergence relation is consistent with 

lim F(Q2, fi^ <oo,gl> 0, A^) = / dp^p{p^ /i^ <oo,gl> 0, A^) = 0. (5.48) 

Q^~*oo JO 

Remarks: 

1. The super convergence relation is the spectral sum rule for the interacting the- 
ory. The superconvergence relation implies the violation of positivity of the spectral 
density of the transverse gluon propagator in Landau gauge. This can be interpreted 
as a manifestation of (gluon) confinement. However, it is not sufficient to conclude 
color confinement. 

2. In 0, it was argued that the superconvergence relation might therefore be 
interpreted as a reincarnation of Haag's theorem; The free theory and the interacting 
theory are inequivalent no matter how small the coupling is. 



6 Infrared fixed point and color confinement 

The color confinement criterion of Kugo and Ojima is concerned with the IR behavior 
of the ghost form factor. Therefore, we focus on the IR behavior of the ghost form 
factor and its relationship with the gluon form factor. 

The gluon propagator for which the superconvergence holds can behave like a 
power of the momenta in the IR asymptotic region. For the ghost, however, the 
superconvergence does not hold in general and hence we can not conclude the ghost 
propagator behaves like a power of the momenta in the IR region. The power-like 
solution of the gluon propagator is possible if and only if the RG function has the 
power series expansion in the coupling constant around an IR fixed point. In what 
follows, we argue that the power like behavior for both propagators is consistent 
with the RG functions which has the power series expansion in the coupling constant 
around an IR fixed point based on the RG equation. 

6.1 IR fixed point 

In what follows, we consider a phase of QCD which is connected to the asymptotic 
free regime, i.e., (3{0) = and P{a) < for sufficiently small a > 0. Let aoo be the 
infrared (IR) value of the running coupling constant «(//), i.e., aoo = Oii^n = 0). 

Suppose that the /3-function is a continuous function of g^. If there is a region in 
which P{a) > 0, there must be at least one such that P{ao) = 0. Therefore, we 
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restrict our consideration to the first negative region of /3(a) < 0, i.e., < a < ao 
where ao is a zero closest to the origin. See Figure |31 If P{a) is negative in the whole 
region a > 0, we can interpret ao = oo. 



P(a) 




uv 



IR. 



V7^ 



a 



a(p) 



a 











JR^^ 








a(p) 







uv 



Figure 5: Two types of the beta function and the corresponding running coupling 
constants. 



Now we proceed to examine which IR asymptotic forms of the propagator are 
possible as solutions of the RG equation. Suppose that the RG functions /3{a),'y{a) 
are single-valued and continuous functions and that they can be expanded into the 
formal power series around a = Ooo'- 

=/3oo + (a - aoo)Pl'^ + " aoofP^"^ + ■ ■ ■ , (6-1) 

7(tt) =7oo + (a - aoo)7(^ + ~ a^^Yli^ + ■ • • , (6-2) 

where we have defined Poo '■= /5(aoo),7oo = 7(aoo)- 

In our argument, the existence of the Landau pole (namely, the running coupling 
a blows up at a certain n = no > before reaching the IR limit /i = 0) is excluded 
from the beginning, since it is considered to be an artifact of the perturbation theory. 
We also excluded the case of Poo = — oo, since the power series expansion around 
a = ttoo loses the meaning in this case. Even in this case, we consider that the 
anomalous dimension 'joo remains finite |7oo| < oo even in the IR limit. Otherwise, 
the field itself will lose its meaning in the IR limit. 

We need to estimate the ratio 7(a)//3(a) to know the IR asymptotic form of 
the propagator where 7(0;) depends on the field in question. We discuss two cases 
separately: (I) Poo = 0, (II) Poo 7^ 0(< 0). In any case, the running of the coupling 
constant is governed by 

Ja{a) P{a) Ja H ' 
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In the case of (I), the integration of ()6.3|) : 

'■"(p) da ^ fP diJL 



shows that the running couphng constant a behaves for small p as 



(6.4) 



(6.5) 



In this case, /3f^ > 0, since /?(0) = 0, /3(aoo) = and /3(a) < for < a < 
Then, the ratio in the case (I) reads 



7(a) 



/3(a) a-ttooPi^ 
which yields 



1 p. 



IR' 



7oo + (a - aoo) ( 7i - ^^ocj^ j + 0{{a - Ooo) ) 



(6.6) 



(6.7) 

after rewriting a{p) in terms of p using ()6.5p . Therefore, the form factor and the 
propagator of the $-field behave respectively like 

F^{p') = {pVaY-, D,ip^):=^^ = {pyaY--\ (6.8) 

In the case (II), the IR behavior of the running coupling constant a(p) is deter- 
mined from 

pip) da ^ rp df,^ 

which is integrated as 

(3^ + jajp) - a^)(3i^ ^(P_Y" ^6 1 

/?oo + (a(a) - aoo)/5f^ \<y) - ^ ' > 

Here the coefficient must be non-positive j3[^ < 0. Otherwise, we run into the 
contradiction. In fact, the RHS goes to zero as ^ for j3{^ > 0, while LHS 
is non- vanishing in this limit a{p) —>■ a^o for j3oo 7^ 0. On the other hand, the ratio 
in the case (II) given by 



7(«) ^ ^ 
/?(«) Po 

yields 



7oo + (a - aoo) (71''' - f^l^'j^j + 0{{a - a^Y) 



(6.11) 



(6.12) 
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This expression converges to a constant in the IR hmit — > 0. Therefore the 
leading IR power behavior of the propagator is still given by the tree form: D^[p^) = 
const. in the case (II). 

Thus we conclude that the second case (II) is excluded for the gluon by the 
superconvergence relation, suggesting the existence of the IR stable fixed point Poo '■= 
/3(q!oo) = 0. Since the /3-function is common to all the fields, the other fields should 
obey the power law ()(i.8j) but with different anomalous dimensions 7^.^^ For the 
gluon, we use 

li = 2K = l (6.13) 

from the above argument based on the superconvergence. The remaining ghost 
anomalous dimension 7^ is determined as follows. 

The non-renormalization [2^] of gluon-ghost-antighost vertex in the Landau gauge, 

i.e., 

Z, = Zl^'Z.Zg = 1 (6.14) 

yields, by taking the logarithm and operating the differential operator /x^, an identity 
among the RG functions 

7» + 27-(a) = ^ = ^<0. (6.15) 
g 2a 



'■oo 



Then we have a relationship in the Landau gauge in the case (I) at a = 

7^ + 27^ = 0. (6.16) 

Therefore, the sign of 7^ is opposite to 7^. Substituting the anomalous dimension 
of the gluon 7^ = 1 into ()6.16|) . we obtain 

7^ = 1, 7g = -l/2 (6.17) 

in the Landau gauge in the case (I). Thus, the gluon propagator is suppressed while 
the ghost propagator is enhanced in the IR region. 

However, this result seems to be incompatible with the superconvergence of the 
ghost in the Landau gauge, which may yield 7^ = 1 for the massive ghost. Note 
that the Landau gauge A = is the fixed point, but it is the IR stable fixed point 
and not the UV stable fixed point. For superconvergence, we need the behavior of 
the UV limit where the perturbation theory is reliable due to asymptotic freedom. 
The superconvergence for the ghost does not hold for A 7^ 0, no matter how small 
A is. The result depends on the initial value of A. Therefore, the treatment of the 

^^Our argument should be compared with that of ^1] (See also Appendix A in for some 
minor corrections) where it is pointed out that the power-series in momentum (without logarithmic 
corrections j^Ol) G{Q^,^^) = ^ndn{g){Q^ / fJ^'^Y'^ can be a solution of the RG equation for the 
ghost formfactor in the Landau gauge A = if and only if the exponents Sn are g independent and 
dn{g) = cxp I ^^ 2[(5n^+7W] y Here the RG functions (3(1) and ^(l) are unspecified in sharp contrast 
with our consideration. The exponent 5 was determined based on the approximate solution of the 
SD equation for gluon and ghost propagators. 
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superconvergence of the ghost propagator is more subtle than the gluon propagator. 
The superconvergence of ghost in the Landau gauge is very special compared with 
that of gluon which has superconvergence irrespective of the choice of A. As we have 
shown in the previous section, this difficulty can be avoided if the ghost does not have 
the massive spectrum and the power-series expansion loses the meaning. In fact, the 
above consideration based on the RG equation suggests that the ghost propagator 
is enhanced in the IR and of the long-range type, while the gluon propagator is 
suppressed in the IR and of the short-range type.^'^ 

6.2 Color confinement 

A sufficient condition for color confinement was proposed by Kugo and Ojima (KO) 
jHHllHn], see also jHOlISI- The criterion in the original form is a condition u{0) = — 1 
for the two-point function of the composite operators: 

^^"^ {a^u - ^) u{p) := J d'xe^^^^-yH0\T{D,C)^{x)g{A, x C)^(y)|0). (6.18) 

In particular, it was shown later that the criterion is considerably simplified in the 
Landau gauge A = jHS]- It was reduced to a condition for the ghost form factor in 
the IR limit: 

lun lQ'AMQ'r' = lun [Gmr' = 1 + m = 0. (6.19) 
Our argument given above leads to the ghost form factor enhanced in the IR: 

G{Q') = {QVa^y- - {QVa'Y'^' . (6.20) 

Therefore, the solution found in this article satisfies the KO color confinement crite- 
rion ()6.19|) . In this sense, color confinement is achieved in the Landau gauge Yang- 
Mills theory. 

As demonstrated in this article, the superconvergence alone is not sufficient to 
draw the conclusion for color confinement. In fact, the KO criterion needs the prop- 
erty of the IR limit ^ of the two-point function for ghost and antighost. In the 
gauges other than the Landau, the detailed information for gluon is also needed to 
verify the KO criterion. According to recent investigations of the coupled SD equa- 
tions for gluon and ghost, contribution from ghost plays the dominant role in the IR 
region. 

This result should be compared with the Nishijima criterion P^94ll95ll^ for color 
confinement 

= 0. (6.21) 

^^As a theoretical possibility, we can not exclude the converse case based on this argument alone, 
namely, the gluon propagator is enhanced in the IR and of the long-range type, while the ghost 
propagator is suppressed in the IR and of the short-range type. However, such a solution with k < 
can not be a consistent solution of the truncated coupled SD equations for the gluon and ghost 
propagators. 
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This criterion is also satisfied in this case, since the combination of the renormahzed 
spectral representation and the renormalized super convergence relation implies 

= / dp^p{p\ /i^ gl, \r) = -, (6.22) 
Jo A* 

as already argued by himself. 

Our argument just given is not complete due to an assumption on the RG func- 
tions, as the following example shows. 



6.3 IR attractor 




Figure 6: The beta function with an IR attractor and the corresponding running 
coupling constant. 



Moreover, there is an exceptional case to the above argument. It is the case of 
IR attractor where (3 function has a pole (not a zero) at a certain point aa, i.e., 
I3{aa — 0) = —00, see Figure IHl In this case, the /9-function blows up (down) at aa 
before reaching the IR limit aoo and hence the /3-function cannot be expanded around 
a = ttoo (Note that aoo > ^a). Some authors [HZj argued that the IR attractor will 
be realized in QCD with small number of flavors, just as the supersymmetric (SUSY) 
gluodynamics has the IR attractor pSjl^ . The exact NSVZ [HE] P function for SUSY 
gluo dynamics is given by 

^^"^ = -2^1-C,(G)a/2.' 

where C2{G) is the quadratic Casimir operator in the adjoint representation, C2{G) = 
Nf. for SU{Nc). Of course, it has the power series expansion around a = 0: 



71=0 



/3(a) = --3C2(G)E(^C2(G)-j . (6.24) 



In the case of IR attractor, the coupling constant a(/i) becomes complex- valued below 
the scale /i^ at which the pole exists = a{fia)- Therefore, we need independent 
study in the case of IR attractor. 
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7 MA gauge 



We have obtained all the possible super convergence relations based on the asymptotic 
behavior of the propagator. See Appendix A for calculations. Summarizing the result, 
the superconvergence relation holds for the diagonal gluon, off-diagonal FP ghost, off- 
diagonal gluon, while the superconvergence does not occur for the diagonal FP ghost. 

In what follows, we choose the gauge parameter j3 = for diagonal components 
and omit to write /5 = hereafter. The the following unsubtracted dispersion rela- 
tions for the propagators and form factors and the superconvergence relations hold 
irrespective of the gauge parameter a for off-diagonal components. The diagonal 
gluons have 

n n 2 2 2 \ f°° J 2 Pdiag {P ^ p j 9 ^ /-r i \ 

Ddiagik ,p ,g ,a) = dp ^2 _ f.2 ' ('^•^) 

17 n 2 2 2 \ f°° J 2P'^Pdiag{p'^-i p'^ig'^iOi) /rjc^x 

Fdiag{k ,p ,g ,a) = - dp ^ , (7.2) 

^—0 P rv 

roo 

0= dk'^pdiag{k'^,P^,g'^,oi). (7.3) 

•J 

The off- diagonal gluons have 

n /, 2 2 2 N p , 2pA(p^/i^c/^a) 

DA{k,p,g,a)= dp , (7.4) 

JO P rv 

p ., 2 2 2 ^ r J 2 P^PA{p^p\g^a) 

Frik ,p ,g ,a) = - dp , (7.5) 



p'^ — k"^ 

0=1 depA{k\p\g\a). (7.6) 



The off-diagonal FP ghosts have 

l\pp[k ,p ,g ,a) = dp — , (7.7) 

Jo P rv 

^ /i,2 2 2 N r J 2 P^PFp{P^P^g^(^) .^gx 

Gppik ,p ,g ,a) = - dp , (7.8) 

J —0 p K 

roo 

0= dk'^pFp{k'^,p^,g^,a). (7.9) 
Jo 

However, the diagonal FP ghost propagator has an unsubtracted dispersion rela- 
tion 



Ad.^,{e,p\g',a)= ^P' ^'^' _p = -^,Gd.ag{k\p\g\a) (7.10) 



,2 ^2 g2 — 1 ^J2Pc{p'^i P^ig'^iC^) _ 1 ^ /j^2 ,.2 „2 

without superconvergence. 

Therefore, the diagonal gluon, off-diagonal ghost and off-diagonal gluon become 
massive. From the viewpoint of Abelian dominance, the mass of the diagonal gluon 
must be smaller than masses m^, mc of the off-diagonal gluon and off-diagonal ghost. 
The off-diagonal propagator must be short-ranged compared to the diagonal one. 
Indeed, a numerical simulation on a lattice indicates this is the case [72]. Moreover, 
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the latest simulation has confirmed this result and provided more detail, i.e., = 
2ma [73] where the gauge parameters adopted correspond to = and Pr = 0. It 
should be remarked that the ghost propagators have not yet been measured in these 
lattice studies [72tl73j . The above results show that the off-diagonal ghost propagator 
We can not say anything about the diagonal ghost. However, only the diagonal ghost 
propagator can have the long-range correlation. This is necessary for the confinement 
to hold, since the only the long-range correlation can convey the information to the 
long distance region where confinement is relevant. 

The massiveness of the diagonal gluon is unexpected at first sight. This is because 
the residual U(l) invariance seemed to impose the masslessness of the diagonal gluon 
and the long-range mode is needed to explain quark confinement, e.g., the linear 
potential. However, we must recall that the ghosts play the role of the longitudinal and 
scalar gluons in the Landau gauge where the longitudinal and scalar gluon disappear. 
Therefore, the linear potential will be provided by the contribution from the diagonal 
ghost. In addition, the massive diagonal gluon after gauge fixing does not contradict 
the general principle of quantum field theory as discussed below form the viewpoint 
of the residual local gauge symmetry. This scenario seems to be not inconsistent 
with the confinement picture obtained in the case of Lorentz gauge fixing. All the 
gluons are expected to become massive, while all the ghosts are expected to be long 
range, provided that the color symmetry is unbroken. In the MA gauge, the color 
symmetry is explicitly broken by the gauge fixing procedure. Therefore, the diagonal 
and off-diagonal components can behave differently. 

8 Conclusion and discussion 

In this article we have examined how the propagators and the form factors for gluon 
and ghost behave as complex-valued functions on the complex plane from the 
viewpoint of general principles (multiplicative renormalizability, analyticity, Poincare 
covariance, spectral condition, completeness) of quantized gauge field theories formu- 
lated in the manifestly Lorentz covariant fashion. In particular, we have emphasized 
that it is the analyticity of the propagator and the corresponding form factor on the 
complex p"^ plane that enables us to connect the UV asymptotic solution into the IR 
asymptotic solution and vice versa. 

We have obtained all the possible superconvergence relations for gluons and ghosts 
in the generalized Lorentz gauge and the modified MA gauge, by making use of their 
UV asymptotic behaviors which are calculable reliably by perturbation theory due 
to asymptotic freedom. The superconvergence relations are obtained by combining 
two sum rules for a propagator and an associated form factor following from the an- 
alyticity on the cut complex p"^ plane. Thereby, we can predict the IR asymptotic 
behavior of the Euclidean field propagator for which the superconvergence relation 
holds. By equating the dispersion relation following from analyticity with the spectral 
representation following from the spectral condition and Poincare covariance, it has 
been shown that the propagator with superconvergence may have massive pole and 
hence become short-range. Moreover, we have discussed the implications of super- 
convergence to color confinement in both gauges. In other words, we have presented 
a scenario of deriving color confinement. 

In the generalized Lorentz gauge, we have confirmed that the gluon has supercon- 
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vergence for arbitrary gauge parameter and hence the gluon propagator may become 
short-range. However, the superconvergence for the ghost should not hold for arbitrary 
gauge. This is consistent with the approximate solutions of the coupled SD equa- 
tions for gluons and ghosts in the Landau gauge []iHllipi2pi6pi7pi8pi9p2np2ip23j . If the 
superconvergence held for the ghost, the ghost propagator would become short-range. 
However, it is reasonable that the IR enhanced ghost propagator be long-range, since 
only the IR enhanced ghost could be a candidate for a carrier of confinement. We 
have concluded that the gluon propagator exhibits the power behavior characterized 
by the IR critical exponent k = 1/2, 

Dt{Q^) = const. + 0{Q^), (8.1) 

while the ghost propagator has a negative and non-integer exponent. 

Supposing the RG functions have the power series in the coupling constant 
around the IR value g'^{fi = 0), we have specified the IR asymptotic behavior for the 
gluon and ghost propagators so as to be consistent with the above consequence of 
superconvergence relations. Moreover, we required that the solution is compatible 
with the existence of the IR fixed point. Then the IR critical exponent of the ghost 
is determined. Thus we conclude that the ghost propagator in the Landau gauge 
behaves like 

AMQ') = {QY'^'- (8.2) 

This result is also consistent with the SD equation approach, if and only if the IR 
critical exponent k takes the value k, = 1/2. The IR asymptotic behavior of the ghost 
propagator ()8.2|) fulfills a sufficient condition of color confinement in the Landau gauge 
presented in j^lj which is originally due to Kugo and Ojima see also jHOlini]. The 
negative half integer exponent for the ghost implies the existence of the branch cut 
beginning at the origin and extending on the positive real axis in the complex k'^ 
plane. This is indeed the situation which is necessary to adopt a new integration 
contour in discussing anew the superconvergence relation for the ghost. 

The IR limit of the running coupling constant g'^{0) can be estimated by solving 
the ghost SD equation by substituting the IR asymptotic solutions ()8.1|) and ()8.2p . 
as already worked out in the SD equation approach, see e.g. Appendix A of ^Tj. In 
Figure [71 we have summarized the relationship between the RG functions (/? function 
and anomalous dimension) and the running coupling constant and propagators. An 
exception to be above argument is the solution corresponding to the IR attractor, 
which has been suggested to exist for QCD with quark fiavors below a certain critical 

^^In the manifestly Lorentz covariant formulation, even the spectral function p{p'^) of the transverse 
gluon is no longer positive. In fact, it is negative for large due to the indefinite metric of state 
space. In view of this, it is difficult to give an intuitive interpretation to the superconvergence 
relation from the physical point of view. In the generalized Lorentz gauge, the longitudinal part 
including the gauge parameter is also renormalized, in sharp contrast with the usual Lorentz gauge. 
Especially, in the Landau gauge, the gluon propagator does not have the longitudinal part. The 
unphysical degrees of freedom such as the longitudinal and the scalar parts are related to the ghost 
and antighost that are also unphysical, but exist even in the Landau gauge. Therefore, we have no 
reason to distinguish the ghost from the other unphysical components of the gluon. 

-'^^If the constant is zero in the gluon propagator H8.1|l . i.e., Dt{Q^) = 0{Q^) [ 0, the the ghost 
propagator behaves like /S.pp{Q'^) ^ {Q^)~^ T instead of H8.2|) . This is nothing but the Gribov 
limit, K = 1 |24j . However, the explicit form obtained by Gribov contradicts with the analyticity 
requirement, since the Gribov form has a complex conjugate pair of poles. 
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number based on the higher order analysis of perturbative series jHZl • This case is to 
be investigated separately. 




Figure 7: The schematic behavior of RG functions, running coupling constant and 
propagators for gluons and ghosts for Yang-Mills theory in the Lorentz-Landau 
gauge, (a) Beta function vs. Running coupling, (b) Anomalous dim. of gluon vs. 
Gluon propagator, (c) Anomalous dim. of ghost vs. FP Ghost propagator. 

Our results do not contradict with the existence of the gluon condensate {s^"^) 
of mass dimension two. The reason is as follows. The superconvergence relation has 
been derived using only the leading term of the ultraviolet asymptotic form of the 
propagator, while the relevant condensate appears as the sub-leading term of the 
ultraviolet asymptotic propagator, as the OPE calculation of the gluon propagator 
shows. Therefore, the existence of such a condensate of mass dimension two does 
not influence the existence or non-existence of the superconvergence relation. [In 
a very special situation (if any) in which all the perturbative terms vanish, such a 
condensate will be directly related to the infrared behavior of the propagator through 
the analyticity.] This is also the case for the gluon-ghost condensate [0310111021 in 
the Maximal Abelian gauge. 

It should be remarked that our results hold for arbitrary choice of gauge fixing. 
Therefore, the IR enhancement of the ghost propagator should hold for both the 
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Landau and non-Landau gauges in the Lorentz type gauge fixing. This is also the 
case for arbitrary gauge fixing parameters in the Curci-Ferrari gauge. Moreover, the 
IR enhancement of the ghost propagator in the Curci-Ferrari gauge is compatible with 
the existence (or finiteness) of the on-shell BRST-invariant gluon-ghost condensate 
■ s^y^ ■\- Xi'^ ■ ^) which is written in terms of the gluon and ghost propagators 
and their form factors as (see ()3.32|1 for the definition) 

\D,,m + AA^p(Q^) 

InQ') + ^FLm + AG(Q2) . 

(8.3) 

For the form factors having the IR power-law behaviors 

F(g2) ^ A ■ iQ'f'^, Fl{Q^) = A' ■ {QY""', GiQ') = B ■ (g')-^ (8.4) 

the integration of the last term is convergent for arbitrary A (after renormalization 
of the UV divergence), if k < 1 (excluding the Gribov limit k = 1 again), since it 
behaves as (e^)^"" for the lower bound e. Therefore, the dominant contribution to 
the gluon-ghost condensate comes from the UV asymptotic region, even if the ghost 
propagator is enhanced in the IR limit. The similar argument can be given also in 
the MA gauge. 

It is also expected that there exists some relationship between the vacuum conden- 
sate and the non-trivial topological object. However, the superconvergence relation 
utilizes only the leading term of the RG improved perturbation theory to character- 
ize the ultraviolet asymptotic behavior which is continued analytically to the infrared 
region. The non-trivial topological object, e.g., instanton clearly appears in the non- 
perturbative treatment. Therefore, it is beyond the present framework of this article 
to extract or separate the effect of the topological object. This issue is to be tackled 
again in the near future. 

In the MA gauge, we have shown that the off-diagonal gluon and the off-diagonal 
ghost have superconvergence and hence all the off-diagonal propagators converge to 
Q^-independent constants in the IR limit. This implies that the off-diagonal modes 
become massive in the MA gauge as is expected from the Abelian dominance 
I7T] . This results do not contradict with the existence of the vacuum condensates of 
mass dimension 2 in the MA gauge inSinZI- We have found that the diagonal gluon 
propagator may he of the massive type, since the superconvergence holds. This is 
a novel feature discovered for the first time in our analytical study, supporting the 
numerical study [73]. This implies that the diagonal gluon becomes massive without 
breaking the residual local U(l) invariance. We wish to emphasize the importance of 
the diagonal ghost propagator in pure Yang-Mills theory with the MA gauge fixing. 
We conjecture that only the diagonal ghost propagator can be short-range to mediate 
the long-range confining force in the manifestly covariant formulation, since only the 
diagonal ghost does not obey superconvergence. These results can be compared with 

^^For the condensate (^^^ • s^^) in the usual Lorentz gauge, the lattice simulation exhibits the 
dominance of the UV contribution for the integral of F{Q^) in the Landau gauge A = 0, see |1()4| . 
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numerical results [721IZS] where the ghost propagators are not measured. There is 
another possibility in the MA gauge, as will be discussed in a forthcoming article. 

In the MA gauge, therefore, the SD equation approach in addition to numerical 
simulations is indispensable to obtain the precise behavior of propagators, especially, 
the diagonal gluon and diagonal ghost propagators as in the Lorentz-Landau gauge. 
In this study, it is very important to reproduce correctly the one-loop perturbation 
result as a self-consistent solution in the asymptotic UV region. Consequently, the 
analyticity allows the UV asymptotic solution to be analytically continued into the 
correct IR asymptotic solution. This requirement will help to constrain the approxi- 
mation (for the vertices) to be adopted for writing down the solvable SD equation in 
MA gauge. At the same time, the truncation of the SD equation must be performed 
in consistent with the multiplicative renormalizability. These points could be good 
criteria of judging the efficiency of the adopted approximations. The SD equations 
of gluodynamics in MA gauge will be discussed in a separate article, see jlHI for a 
preliminary result. 

From our viewpoint, recent studies of the coupled SD equations for gluon and 
ghost propagators in the Landau gauge performed by Bloch [TBlll7ll23j are very inter- 
esting in the following respects. He has proposed the truncation scheme respecting 
the multiplicative renormalizability and reproducing the leading order resumed per- 
turbative results (beyond one-loop result) in the UV region ^H]- Then he has actually 
solved the coupled SD equations for gluon and ghost propagators by including for the 
first time all the diagrams (including two-loop diagrams in the gluon equation) and 
has found the numerical solution over the whole range of momenta characterized by 
the power behavior in the IR region |2S]- Among other things, his solution includes 
the case of k = 0.5, although all the approximate solutions obtained so far has k 
slightly larger than 0.5. His solution seems to be an approximate solution which is 
consistent with the general principles of gauge field theories. 

Here we comment on whether or not our results are affected by the Gribov prob- 
lem. Gribov j21j studied how the requirement of preserving positivity of the ghost 
propagator influence the original Yang-Mills theory. This requirement forces one to 
restrict the region of functional integration over the gauge fleld to a subspace con- 
nected to the origin (called the flrst Gribov region) inside the Gribov horizon. Our 
results still hold even in the presence of Gribov problem, as far as the analyticity is 
kept even after restricting the functional integration. Otherwise, it becomes impos- 
sible to perform the analytic continuation of the Euclidean propagator back to the 
Minkowski region. 

In order to know the analytic structure of the propagators more explicitly, it is 
desirable to solve the SD equations on the complex plane. Such investigations 
have been performed only for the quark SD equation |100[I101] . Solving the cou- 
pled SD equations for gluon and ghost on complex plane will be a challenge in 
future. See jl()2j for numerical simulations. In a subsequent article we will discuss 
the effect of including the fermion flavors and clarify the interplay between color 
confinement/deconfinement and chiral symmetry breaking/restoration. 
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A Calculation of UV asymptotic behavior of the 
propagators 



A.l Generalized Lorentz gauge (Curci— Ferrari gauge) 

The running of the gauge coupling constant is governed by the differential equation: 



9/i 
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and the gauge fixing parameter are governed by the differential equations: 
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Eq. ()A.1|) is a closed equation and it is solved exactly. 



22 C2(G) 1 M 
3 (47r)2 "^Aqcd 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



1 + ftf^lMln^^ 

where we have used the boundary condition gn^^o) = oo at /io = Aqcd- The remain- 
ing two equations ()A.2|) and ()A.3|) are coupled equation. Using the solution ()A.4|) and 
rewriting the derivative -^jJ^-i- in (IA.2|) . (1A.3|) as 
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we can eliminate dependence in flA.2|l and ()A.3j) as 
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They are still coupled equation and we can not solve them exactly. Therefore, we 
restrict our consideration to the subspace in which the flow is always confined. It 
is the fixed subspace consisting of the fixed points of the RG. The fixed point (to 
one-loop order) is obtained by solving the algebraic equation simultaneously: 
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Figure 8: RG flow in the A) plane. 
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Figure 9: RG flow in the (a, a') plane. 



We find one fixed line consisting of fixed points and three isolated fixed points in 
the (^R, Ar) plane, or equivalently four isolated fixed points in the (or, a'^) plane: 
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A. The line of fixed points: Ar = 0, .^r G R corresponds to an isolated fixed point 
(aR,0 = (0,0). 

B. (^R, Ar) = (|, f ) corresponds to (or, a'^) = (f , 0). 

C. A) = (0, y) corresponds to (or, c^r) = (0, y ). 

D. (^R, Ar) = (1, y) corresponds to (or, Or) = (y , -y) 

If the two parameters .^r, Ar are set equal to one of the fixed points, the theory 
remains forever on the fixed subspace. If the system starts from other points and the 
scale /i is decreased, it evolves into the infrared (IR) region according to a couple of 
differential equations ()A.1|) . ()A.2|) and ()A.3|) . 

Among the RG flows, the five RG flows (I, E, IE, IV, V) connecting the fixed points 
A,B,C,D form the watershed (or backbone) in the flow diagram. 
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(A.9b) 
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Figure 10: The flow of the gauge parameter A in the generalized Lorentz gauge. The 
arrow is directed to the IR region, i.e., decreasing the renormalization scale /x. 



In what follows, we use A to denote the initial value, X := Xr the running parame- 
ter and A* the fixed point of RG. First, we consider the fixed subspace II, i.e., ^r = 



in which A* := y and the Ar is governed by 

^Ar _ 3 
ainlnx^ 22 

This equation is easily integrated 
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Note that A = for A = and hence A = is an IR fixed point, while the point 
A = A* is an UV fixed point. A = A* for A = A*. For other initial values of A, we find 
that A I A,, for A > A* and A t A* for A* > A > 0, while A — ^ — oo for A < 0. 
The anomalous dimension of the gluon field is ^-independent and is given by 



1a{91) = -1{X.-Xk)C,{G)-^^, 

and the ghost field has the (also .^-independent) anomalous dimension 
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Substituting ()A.13|) for A = Ar into ()A.14|) and ()A.15|) . we can regard the anomalous 
dimensions as the functions of the gauge coupling g"^ alone. 
For the A in the neighborhood of A*, ()A.13|) reads 
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Therefore, the exponential factor for the gluon reads 
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where we have used 1 + 4^ = 1 + \ = j-. For the ghost, the factor reads 




The first term diverges as g ^ 0. In this case the superconvergence relation does not 
exist. In this case, the asymptotic form of the ghost form factor is given by 

G{e) ^ ^1 + (\n ^oo as ^ oo. (A.19) 

The ghost spectral function has the asymptotic form 

P.P(P). (l + l^) (in g*". (A.20) 

These results agree with those of j56j . 

For sufficiently small A in the neighborhood of A = 0, ()A.13|1 reads 

_13 

A = A(^j (A.21) 
Therefore, the exponential factor for the gluon reads 




This result is the same as the above result ()A.17j) . For the ghost, however, the factor 
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reads 




This factor diverges for non-zero A as (7 ^ no matter how small A is chosen. In this 
case the superconvergence relation does not exist. Only for an initial value A = 0, 
this factor goes to zero as (7 — *■ and the superconvergence relation holds. Therefore, 
there is a discontinuity at A = for the ghost factor. For the superconvergence to 
be realized for the ghost propagator, we need a fine tuning of A. In this case, the 
asymptotic form of the ghost form factor is given by 

In^ ^00 as fc^^oo. (A.24) 
The ghost spectral function has the asymptotic form 

——1 

pMe) = -{kr'^^{\n^Y . (A.25) 

This case has not been worked out in 

Next, in the fixed subspace I, i.e., ^ = 1/2, the final result is obtained by replacing 
II by II with the identification A* = y. 

13 \ -1 /■ , 13 \ -1 

The results for the A in the neighborhood of A* are obtained by replacing the exponent 
^ by The results for sufficiently small A in the neighborhood of A = are obtained 
by replacing the exponent || by j^. Finally, in the fixed subspace III, i.e., ^ = 1, the 
same result as I is obtained, since the theory has the refiection symmetry ,^ ^ 1 — 
They are new results for the generalized Landau gauge. 

Thus the ghost can escape from the superconvergence except for the Landau gauge 
A = 0. This point is very important for our result due to general principle to match the 
results of approximate solutions of the Schwinger-Dyson equation for gluodynamics. 
The Landau gauge limit A | may be rather subtle. 
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A. 2 Modified MA gauge 

We discuss only the case G = SU{2). It was shown giEHEniEHl that the RG 
equation for the gauge fixing parameter /3 of the diagonal part is given by [36tl37] 
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and the gauge fixing parameter a of the off-diagonal part obeys the RG equation 
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We apply ()A.5|1 to rewrite the differential equation (jA.27j) into the form which 
does not explicitly depend on g'^: 
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which is easily solved: 
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Figure 11: The flow of the gauge parameter a in the modifled MA gauge. The arrow 
is directed to the IR region, i.e., decreasing the renormalization scale 



In the similar way, ()A.28jl is cast into 
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First, we consider sufficiently small a, neglecting the order term, 
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which yields 
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Next, we take into account the O(a^) term too. Applying the formula 
dx 2 2ax + b 
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This shows that a —oo as /i ^ oo irrespective of the value of a. Note that arctan a; 
has multivalued, unless — 7r/2 < arctan x < 7i/2. 

The anomalous dimension of the diagonal gluon is given by |H7j 
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while the off-diagonal gluons have the anomalous dimension |^IH7] 
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The anomalous dimensions of the diagonal ghost and diagonal antighost is calculated 
as I3II 
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while the off-diagonal ghosts read 
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For the derivation of anomalous dimensions, see P7| l36j. 
For the diagonal gluon, the correction factor reads 
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The super convergence relation holds for the diagonal gluon. 
For the off-diagonal gluon, 
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for p <0 (as ^ ^ 0), 
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where we have used 
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The RHS oscillates, but goes to zero for arbitrary a for /5 < 0, while it goes to infinity 
for j3 > 0. The super convergence relation holds for (3 < 0. 
For the off-diagonal ghost. 
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As fi ^ oo {g ^ 0), the RHS goes to zero for arbitrary a for /3 < 0, while it goes to 
infinity for j3 > 0. The super convergence relation holds for (3 < 0. 

For the diagonal ghost, the anomalous dimension of the diagonal ghost has the 
same magnitude as that of the diagonal antighost with opposite sign, i.e., 7c(5'r) = 
— 7c((7|.) and hence the corrections cancel. 




There is no correction to the tree form for the diagonal ghost-antighost propagator. 
Therefore, the superconvergence does not occur for diagonal ghost and antighost. 
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